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Abstract 


A  data  bank  system  based  on  the  concept  of  data  bins  is  described 
and  examples  of  its  use  with  gravity  and  topography  data  are  given.  This 
concept  allows  for  the  efficient  storage,  retrieval  and  manipulation  of 
large  amounts  of  geophysical  data  and  is  applicable  to  other  data  sets, 
for  example,  heat  flow,  magnetic,  isopach,  and  geochemical  sampling 
data. 


Land  gravity  and  topography  data  currently  held  at  Lamont  is  being 
reformated  into  data  bins  defined  by  geographical  1°  x  1°  squares.  With 
time,  bathymetry  and  marine  gravity  data  will  also  be  reformated  into 
these  bins  creating  a  global  gravity/ topography  data  bank. 

Each  sampling  position  within  a  bin  is  stored  in  card  image  format, 
one  card  image  per  station  (or  data  point),  and  includes  the  coordinates  of 
the  station,  height  above  sea  level  (or  depth),  observed  gravity,  free- 
air  gravity,  bouguer  gravity  and  source  information.  If  other  geophysical 
data  exist  at  this  locality,  for  example,  heat  flow,  the  record  is 
easily  expanded  to  include  this  extra  information. 


Introduction 


Plate  tectonics  is  a  unifying  working  hypothesis  which  provides  a 
kinematic  model  of  the  upper  layer  of  the  earth  and  a  global  framework 
for  geologic  processes.  It  can  be  used  to  make  qualitative  predictions 
about  most  phenomena  studied  by  the  different  disciplines  of  the  Earth 
sciences.  This  leads  to  a  re-examination  of  large  amounts  of  geological, 
geophysical  and  geochemical  data  already  acquired  in  order  to  formulate 
precise  theories  for  earth  processes. 

The  plate  tectonic  hypothesis  addresses  problems  on  a  global  scale 
and  so  to  integrate  existing  data  for  use  in  a  global  sense,  major  restruc¬ 
turing  and  reformating  of  this  data  is  usually  necessary.  Most  data 
acquisition,  storage,  and  retrieval  systems  have  generally  been  inherited 
from  the  exploratory  investigations  of  ocean  basins  over  the  last  decade 
in  that  they  utilized  time  sequential  (or  cumulative  distance)  format. 

While  this  is  acceptable  for  low  density  sampling  and  storage,  this 
format  rapidly  becomes  cumbersome  for  higher  density  sampling.  Meaningful 
data  retrieval  is  a  need  which  is  becoming  more  evident  as  our  ability 
to  gather  and  store  data  becomes  increasingly  sophistocated . 

What  follows  is  the  processing  approach  taken  by  the  Lamont  Gravity 
Group  in  utilizing  the  pre-existing  gravity  and  altimetry  data  banks  in 
global  gravity  studies. 

Data  bank 

The  acquisition  of  gravimetric  and  topographic  data  by  various 
governmental  and  scientific  institutions  over  the  last  decade  has  provided 
the  basis  for  a  global  gravimetric  and  topographic  data  bank  (once  the 
data  has  been  made  compatible  in  both  the  method  of  reduction  and  the 
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storage  format) . 

As  the  sampling  of  data  approaches  that  of  an  areal  coverage 
(rather  than  line  coverage),  it  is  reasonable  to  sub-divide  geographical 
areas  into  a  number  of  compartments  or  bins,  the  boundaries  of  which 
form  an  integer  grid  of  latitude  and  longitude.  With  the  world's 
surface  binned,  any  gravity  station  will  naturally  belong  to  one  bin. 

The  purpose  of  the  bins  is  to  reduce  the  amount  of  data  to  more  manage¬ 
able  working  sizes  such  that  searching,  sorting  or  updating  of  the 
data  can  be  done  efficiently  and  subsequent  processing  more  direct. 

Any  size  data  bin  can  be  created,  for  example,  one  degree  by  one 
degree,  ten  degree  by  ten  degree  etc.,  but  the  size  selected  will  be 
dictated  by  the  sampling  density  of  the  data  over  large  areas.  Because 
of  the  detailed  coverage  of  land  gravimetric  data,  we  have  decided  on 
one  degree  X  one  degree  bins.  Again,  the  selection  of  bin  size  is  a 
comprimise  between  the  number  of  points  within  any  one  bin  (which  control 
the  efficiency  of  searching/sorting)  to  the  total  number  of  empty  bins 
created  for  an  area.  The  bins  are  stored  on  computer  disk  with  a  unique 
filename  created  from  the  latitude  and  longitude  of  the  bin's  bottom  left 
hand  corner  (both  northern  and  southern  hemispheres).  The  naming  conven¬ 
tion  is; 

[B]  [H]  [LLNNN]  [IDEG]  .  [NIB] 

where, 

B  -  bin;  if  more  than  one  data  set  is  to  co-exist  on  the  same 
mass  storage  device,  this  letter  could  indicate  the  type  of 
data, 

H  -  hemisphere,  N  for  north,  S  for  south, 

LLNNN  -  LL  is  the  latitude,  NNN  is  the  longitude  of  the  bottom 
left  hand  corner  of  the  bin. 
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I  DEG  -  the  bin  size  in  integer  degrees, 

.NIB  -  PDP  file  extension. 

In  creating  the  data  bins,  it  has  been  assumed  that  latitude  is 
defined  from  -90°  to  +90°,  positive  being  in  the  northern  hemisphere 
(0°  being  positive),  and  longitude  is  defined  from  000°  to  360°. 

(negative  longitudes  can  be  used  as  input  but  are  internally  converted 
to  the  equivalent  easting).  The  coordinates  of  any  gravity  station 
falling  within  the  bounds  of  a  defined  bin  are  appended  to  the  file 
along  with  other  pertinent  information  in  card  image,  PDP  sequential 
file  format.  By  using  these  files,  the  PDP  system  routines,  such  as 
editing,  printing  etc.,  are  still  available  to  the  user. 

The  card  image  contains; 

Latitude,  longitude,  observed  gravity  (absolute,  mgals), 
topography  (metres,  +ve  ASL),  free-air  gravity  computed 
by  using  the  1967  International  Gravity  Formula,  bouguer 
gravity,  various  processing  flags  and  finally,  alphanumeric 
information  relating  to  the  source  of  the  data. 

The  purpose  of  the  source  data  is  to  include  an  option  for  selective  pro¬ 
cessing  of  the  binned  data.  For  example,  data  can  be  preferential ly 
selected  or  rejected  depending  on  the  source  data  during  the  contouring  of 
the  data.  Because  of  the  absolute  size  and  number  of  files  contained  in 
the  data  bank,  it  is  stored  on  a  dedicated,  demountable  disk. 

Programs 

Two  main  programs  presently  exist  for  creating/updating  the  data 
bank  and  manipulating  the  bank  to  create  files  compatible  with  pre-existing 
Lamont  programs.  Program  BINNER  creates  and/or  updates  the  data  bank. 
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Program  PBGT  is  an  example  of  the  way  in  which  the  binned  data  can  be 
automatically  accessed  -  in  this  case,  to  create  profiles  of  topographic 
and  gravimetric  data.  Other  programs  that  will  be  added  to  the  system 
in  time  are; 

ABGT  -  annotate  bathymetric  and  gravimetric  values  in  a  given  area, 
area, 

CBGT  -  automatic  contouring  of  bathymetric  and/or  gravimetric 
data  in  a  given  area, 

IBGT  -  Interpolate/average  data  within  bins, 

DBGT  -  Dot  or  station  plot  of  bathymetric  and  gravimetric 
values  (for  example,  figure  1). 

(i)  BINNER 

Program  BINNER  bins  raw  gravity  and  topography  data  once  the  data 
has  been  pre-processed  to  compatible  formats  etc.,  from  external  tapes. 
Examples  of  its  use  are  shown  on  pages  B1  -  B3  for  geographical  coordin¬ 
ates  which  span  both  the  equator  and  the  Greenwich  meridian.  The  various 
responses  are  more  fully  defined  in  the  Input/Output  section.  The  real 
value  of  binning  the  data  becomes  apparent  when  it  is  necessary  to 
locate  data  in  various  geographical  areas,  for  example,  in  producing 
contour  maps  or  constructing  data  profiles.  Once  the  data  is  binned, 
subsequent  routines  must  be  able  to  automatically  open  the  appropriate 
bins  to  access  the  data. 

(ii)  PBGT 

Say  that  the  data  in  a  particular  geographical  area  has  already 
been  binned  and  we  wish  to  project  the  data  onto  some  baseline.  We 
define  an  azimuth  by  the  initial  and  final  coordinates  of  the  baseline 
and  a  distance  perpendicular  to  the  baseline  for  which  data  can  be  used 
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. 
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in  the  projection.  We  could  sequentially  search  through  the  bins  testing 
whether  each  point  is  contained  within  this  "window"  -  a  very  time 
consuming  operation  (as  well  as  defeating  the  purpose  of  binning). 

Instead,  program  PBGT  calculates  a  subset  of  bins  which  may  be  used  and 
then  opens  those  bins  through  which  the  baseline  passes. 

The  rectangle  defined  by  the  initial  and  final  coordinates  of  the  base¬ 
lines  defines  the  subset  of  bins  that  may  potentially  be  used  in  constructing 
the  profile.  The  X-Y  coordinates  for  each  bin's  vertex  relative  to  the 
initial  coordinate  of  the  baseline  is  calculated  and  if  there  is  a  change  in 
sign  in  the  Y  value,  the  baseline  must  pass  through  this  bin  and  so  it  is 
opened  and  sequential ly  searched  for  points  contained  within  the  profile 
wi ndow. 

The  output  from  this  program  is  in  Lamont  project  file  format  (see 
Brown  book  program  PRJCT)  so  that  it  can  be  readily  used  by  other  existing 
Lamont  programs. 
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Input/Output  control. 

PROGRAM  BINNER 

(1)  Input  geographical  data  limits  remembering  that  latitude  is  defined 
as; 

0  <:  lat  £  90  Northern  hemisphere 

-90  <:  lat  <  0  Southern  hemisphere 

longitude; 

000  ^  Ion  ^  360. 

The  geographical  coordinates  WEST,  EAST,  SOUTH,  NORTH  are 
equivalent  to  LEFT,  RIGHT,  BOTTOM,  TOP  of  a  map. 

(2)  Input  the  size,  in  integer  degrees,  of  the  bins  to  be  produced, 
i.e.  1  x  1°,  10°  x  10°.  For  most  geophysical  purposes,  1°  x  1° 
and  10°  x  10°  (100  x  100  km  and  1000  x  1000  km  respectively) 
should  adequately  represent  data  distribution  density.  There  is 
obviously  a  compromise  between  the  number  of  bins  stored  on  disc 

to  the  maximum  size  any  one  bin  can  become  before  it  is  inefficient 
to  search  the  bin. 

(3)  Input  processing  flag;  if  the  bins  already  exist  and  new  data  is 
being  binned,  answer  is  0  (zero  for  no).  However,  if  a  new  area 
is  being  binned,  answer  is  1  (one  for  yes). 

The  program  now  responds  by  listing  the  input  data.  Note  that  the  latitude 
and  longitude  differences  are  incremented  by  one  block  ensuring  that  there 
will  be  sufficient  bins  for  the  end  points  (for  example,  binning  from 
(+10,  145)  to  (+25,  180),  the  last  point  (+25,  180)  belongs  to  the  bin 
BN25180,  not  BN24179  and  may  well  be  the  only  point  in  this  bin). 

Each  bin  name  is  displayed  on  the  console,  if  the  file  already  exists, 
the  original  bin  remains  unaltered,  the  program  continuing  to  the  next  bin. 
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(This  implies  that  there  are  no  other  file  versions  besides  . N I B ; 1 ) .  If 
the  file  does  not  exist,  a  file  with  that  bin  name  is  created.  After  all 
the  bins  have  been  created  for  a  geographic  area,  the  raw  data  file  name  is 
entered  so  that  binning  can  begin.  Each  point  is  read  in  turn,  and  from 
its  coordinates,  the  corresponding  bin  is  opened  and  the  record  is  appended 
to  the  file.  This  bin  remains  open  and  data  is  appended  until  a  station 
outside  of  the  bin  is  encountered.  Note,  there  is  no  check  on  duplicate 
data  points. 

If  a  raw  data  point  is  read  for  which  a  bin  does  not  exist,  the  entire 
card  image  is  written  to  a  file  called  BINMISS.DAT.  Sometimes  the  raw  data 
file  may  contain  various  areas  and  only  a  portion  needs  to  be  binned  in 
which  case  the  BINMISS  file  is  redundant.  However,  if  the  missing  bins 
resulted  from  an  oversight  by  the  user,  the  BINMISS  file  can  be  input  to  the 
binning  routine  to  re-bin  the  skipped  data  once  the  missing  bins  have  been 
created. 

As  an  example  of  the  entire  binning  process,  consider  the  sampling 
distribution  of  gravity  an  topography  stations  shown  in  figure  1.  One  degree 
x  one  degree  bins  would  subdivide  the  data  as  indicated.  The  raw  data  file 
for  this  area  was  called  USHRT.FAE  and  was  stored  on  DB2.  Program  BINNER 
was  used  to  bin  this  raw  data  as  shown  on  page  B1 .  To  check  the  bins,  the 
system  DIRECTORY  instruction  was  used  as  shown  on  page  B5.  Notice  the 
empty  bins  (e.g.  BN4923801 .NIB)  and  their  position  in  the  geographical  grid  - 
they  represent  bins  for  points  that  may  have  a  latitude  of  49°  or  a  longitude 
of  238°.  Now  that  the  data  is  binned,  we  wish  to  construct  a  profile  of 
the  gravity  and  topography  data  from  (48,235)  to  (49,242)  with  a  20  km 
window.  The  resultant  projected  files  (with  PRJCT  format)  have  been  plotted 
on  the  console  and  are  shown  in  figure  2. 
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LN3  =  .4Am£2-(  L.4 1  *10o>LN2*  1  0) 

ENCODE  C 15.90 , IFILE)  HEM( I HEM) , LT 1 , LT2 , LN 1 , LN2 , LN 3 , 1  DEG  1 , 1 DEG2 
90  FOR*  AT  (  'd',Al,7Il,'.Nld;i') 

oRITECS,  10 0)  JC0U.4T,  (  IFILE(K)  ,  K  =  1 ,  15) 

100  FORMAT ( 2X , 2 H (  ,I4,2H  ),15Ai) 

I  FILE ( 30 ) aO 

OPEN ( UNI 1=1 , NAMES IFILE, TYPE= 'OLD'  , ERR  =  110) 

GO  TO  120 

110  OPENCU  UT  =  1  ,NAME  =  I  FILE,  TYPE= 'NEW  ,  ERRs  120) 

120  CLOSE ( U N I T= l ) 

C 

C  READ  RAa  INPUT  DATA  AND  elN  INTO  AbGVE  FILES. 

130  *R ITE l 5 . 140) 

140  FORMAT!'  TYPE  RAN  DATA  FILE  NAME  -  30A1  ') 

ACCEPT  150,  IFILE 
150  FORMATC30A1) 

IFILE(30)=D 

OPEN ( UNITs2, NAME= IF LLE, TYPES' OLD'  , ACCESSs • SEQUENTIAL ' .ERR  =  130) 
OPEN  1  UNITS  3 , NAMES' SINMISS.DAT '  , TYPE= *  UNKNOWN '  ,  ACCESSs » APPEND '  ) 

C 

loO  READ! 2, 170,EN0=200)  RL A T , KLON , GOBS , TOP . FA  A , BOUG , ( ITC C I ) ,1=1,8), 
1  (SOUkCEC 11) , 11=1 , 20) 

170  FURMAT(2F10.4,2F10.1 ,2F6.1 ,811,20 A  1 ) 

IHEMsl 

IFCRLAT.LT. U.O)  IHEMS2 
IFCRLON.LT. 0.0)  RLON=360 . O+RLON 
C 

C  uPEN  FILE  RELATING  TO  READ  DATA  POINT. 

LAT  sABS(HLAT) 

L 0  N  =R  LON 

IFC AMODCRLAT,DEG)  .EQ.O.O)  GO  TO  175 
IFCRLAT.LT. 0.0)  LAf=LAT+lUEG 
175  IFC1HEM.NE.IHM0LD)  GO  TO  180 
IF  t  LAT  . N£. LA  TOLD )  GU  TO  180 
IFC LON  .EQ.LONOLD)  GO  TO  190 
180  CLOSE  C  U M T= 1 ) 

C 

NEARSsCLAT/CDEG)*IDEG 
NEAR  n=(LU  14/10  EG  )*IDEG 
C 

LT 1  = NEARS/10 

LT  2  sNEARS-LTlMO 

LN 1  sNEARX/100 

LN2  SNEARW/10-L  4  1*10 

LN  3  sNEAR*-CLNi*100  +  LN2*10) 

ENCOOEC 15,90, IFILE)  HEM U HEM) ,LTl  , LT2 , LN 1 , LN2 , LN 3 , 1DEG 1 , IDEG2 
IFILEC30)=0 


■ 


> 


- 
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KOkl'KAiN  TV-PoUS  «/02-SlE 
dlNutk  .  FTN  /  Tk  :  9lOCnS/*K 


23-GC  r-HU 


0  U  9  S 
00k  o 
0047 
OuSH 
ooyk 

0100 
010  1 

0102 

0103 

0104 

0105 

OlOb 

0107 
0  1 0  3 


OPEN IUN 1 F=l , n4*E=IE1LE,) 1PE= 'OLD ' , ACCESSs ' APPEND • ,KRR=300 ) 
140  LaIuLD  s.NEARS 
Li  J  n  '1 LO  =  i  E  A  R  •» 

lHNULl)  =  1  iER 

*KI  rt(  1 ,  1  70)  RLA  r,ku.0N  ,Gu»S,TCP,EAA,6aUG,  (ITCCi)  ,1  =  1  ,8)  , 

1  (iOURCEili) ,ii=l ,20J 

GO  TO  160 
200  CONTINUE 

C(jU3E(  UNIT=l  ) 

CL.03EC  UNIT  =  2  ) 

CLOSE ( UN  IT* 3 ) 

CALL  EXIT 

300  *RlTE( 3, 170 )  RL A T , k LOn  f GObS , TOP , E A A , BQUG , ( I TC ( I ) ,1  =  1,8), 

1  (SJUKCECli) ,11=1 ,20) 

GO  TO  IbO 
ENU 


A5 


PROGRAM  PBGT 


r  up  I'hMfj 
*-ovjX.It’T« 


u  o  u  2 

UOilJ 

ooo4 
ooOb 
OOob 
ooo  l 
UV  U8 

u»u4 

oo  10 


uu  1  i 

Ov  1  2 

oo  1  j 
ou  m 

Ov>  15 
00  10 


o><i7 

OO  lo 

ooi4 
oo  2U 
00*1 
uu22 
002  3 
oo  2*t 
oo2b 
oo^o 


oo27 

oo2» 

0024 

oo  jo 

00  J1 

OO  J2 

00  jj 
oo  j4 
oo  jb 
oo3t> 
oo  j  / 

OO  jb 

00  jy 
0040 

0041 

0042 


0043 

OO44 

004b 

0  O  4o 

004  1 

0040 

0044 

OobO 


IV-Pf.Us  v  <•  2-b  1  K  1  4  :  2  1  :  b  u  28-JAW-81 

/  T H  :  PLUCKS/ WR 


PACK  1 


(J»4 

C*  + 
C** 

c** 

c** 

c** 

c** 


c 

c 


c 

c 


c 

c 


PRuGrmM  P6GT 

(P)RUF1L£  (U)l  N\EO  (G)RAVITY  AmD  C  T  )  OPCGR  A  PH  1 C  DATA. 
ROUTINE  13  t'0  0  1  v  ALiEN  f  iu  PROGRAM  GRAPE.  lNPu'1'  TO  THIS 
PROGRAM  IS  HINNED  DAT*  PREVIOUSLY  COnSTKUCIFU 
USInG  PROGRAM  ti 1 N w t H  . 


rtUThOR 


G  .  D  .  X  rt  R  N  E  A  (APKIL,1980) 


10 

2o 

3  0 


40 

bo 
0  0 


LYTE  LAdELl (30  )  ,LAdEL2(20) ,LABEL3(20) 
hYTE  H E (  2  ) 

u  lMr>SIUN  PUS  (  2000)  ,  G  R  V  (2000)  ,  l'UP(20QU) 

I  oTtGEK  DIVIDE 
REAL  N  I f! E 

DATA  NPTS,LU,<1,LUN2,RP0INT,MNE,HEM/0.1 , 2  . -0 . 1 . 0 . 99  999  ,  '  N  '  .  *  S  '  / 
uaTA  D  EG  K  a  D / 0 . 0 1 7  4  b  3  2  9  2  b / 

CALL  EK RSET( 29, .TWO E. , .FALSE. , .TRUE. , .TRUE. ) 

CaLo  ERKSEi'C  3 0,. TRUE.,. FALSE.,. ThUE.,. TRUE.) 

REAL  JUTPUT  DATA  FILE  NAMES . 

UKirElS,2o) 

kokmatc  enter  uutpur  filenames,  2oai,2oai  ') 

Rt'AU(b,3o)  LA3EL2 ,  LABEL3 
FUR  1  <* T  (  2  0 A  1  ) 

LABEL2 ( 20 ) so 
L A  BtL3 ( 2  0 ) so 

READ  CUURDlN  AXES  OF  ST  ART/END  OF  PROJECTION  LINE. 

«R1TE (5,40) 

FUR  MAH'  ENTER  PROFILE  COORDS  -  WL  AT  ,  WLOn  ,  FLAT  .ELCn  ') 

ACCEPT  * , n L  A  T , *  L  J  N , E  L  A  T , E  L  U  N 
.«HITE(  5, 50  ) 

FiiR  4  a  r  (  •  enter  projection  origin  -  klat.hlom  ') 

ACCEPT  *,KLAT, RLON 
«RITE(b,t>0) 

FORMAK'  EJTtR  aInuUW  LENGTH  (KMS)  ') 

ACCEPT  *,  aINDUw 
4 1  n  Du  w  sW  I N  DO .« /  2 . 0 

CALCULATE  AZIMUTH  AND  LENGTH  OF  THE  PROJECTION  LINE. 

CALL  AZUI5TULAT, «LON , ELAT , ELUN , R A NGE , A ZMUTH ) 


I  HE. msi 

IF ( nLAT.LT. 0.0) 
IFULON.LT. 0.0  ) 


I H  E  M  s  2 
WL0W=J60.U+wL0N 


bb 


IF  (ELON.LT. 0.0)  ELOns3qO  .(,)•*■£  LON 
IF (KLuN.lI .0.0)  RLON=JbO .O+RLON 


Li  1 N  m s  a  LON 
LUNE=ELGN+N INF 
L  A  T  N'  =  A  r)  S  (  >•  L  A  T  ) 
L AT£=ABS ( ELAT ) 
'i  a  1 T  E  (  5,6b  ) 
FORM  A  T  l  1  N  X  N 
ACCEPT  *,  IDEG 
DEG  =IDEG 


HI, OCR  SIZE  (INTEGER  DEGREES)  ') 


7  u 


>,  KITE  (5, 70)  w  L  AT,  w  LON, ELAT, ELON,RANGE,AZMUTH,  WINDOW,  LA  BEL2,LAbEL3 
FQRMhT(/,SX, 'PROGRAM  PBGT  -  PRODUCTION  OF  PROJECTED  LAND  DATA',//, 


1  IX, ' PROJECTION  LINE  COORDINATES 

2  IX ,' PROJECTION  uENGTH  (AMS) 


' ,4(1X.F7.2. IH. ) ./, 
•  ,3X,F8.2,/, 


3  IX, 'LINE  AZIMUTH  AND  HALF  WINDOW  LENGTH  -  ' , 2 ( 2 X , F6 . 2 , 1 H , ) , / , 


4  IX, 'OUTPUT  DATA  FILES  LABELLED 

IF ( A«Uf)  (  WLAT.DEG)  .  EQ  .  0  .0  )  GO  TO  80 
I  FULAT.LT. 0.0)  LA  TWsLAT  «  +  IDEG 
80  I r  ( A nUO ( ELAT , u£G ) .EQ.O.U)  GO  TO  90 
IF  (  EL  AT.  LT.  0.0)  L  A  TE  =  l,  ATE  +  I  DEG 
LUnW  s 1 LOn * / I  DEG ) * I  DEG 
LA  Tw  s(LATw/IO EG )*IOEG 

90  ISIGN  sSIGW(1.05,*LAT) 

IS  IGN2=SIC,N(  1 .05,  EL  AT) 


-  '  .20A1  ,/,39X,20Al  ,/) 


.",7 


' 


- 


■ 


■ 

■ 
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uusi 

v  b  2 

uv>5i 

v>Ob4 

o  Ub5 

UObo 


Obb7 

VW5# 

o«bi» 

l/Ubl; 

Ul70l 


v  v  o  2 
v  ^  b  J 


uw03 

o  <y  o  O 

U  o  o  / 
'-'Woo 

lyu  /o 

v  v  /  1 


w;  / 1 
uu  7  j 

o  u  /  4 

vi«  7b 


u  U  /  o 
u  0/7 
v>o7  b 

m  /  y 

O  O  0  tj 
U0b  1 

ooo2 

OOob 

Oo  o4 

oo<3  b 

UUBO 

o  o  o  7 
ooao 

Ov>0  ^ 

uo^o 

oo^i 
ooy  i; 
uo^j 


uo^4 


0^55 

ooy0 


0  0  y  7 
00^8 


LATul  =LAr**iSiGN 

Li  A  i  U  =A  AS  (  L.  A  T  m  *  IS  I GO- LA  IE*  ISI  GN  2  )  +  I  DEG 
LuOD  =  ABS  (  L« i J 1 4  W  —  { » f .) ■  < tl  )+IDEG 

I F  (  uUN  a  .  GT  .  i<t>  0  &  )  bO  N  U  =  3  b  0  —  LUND+IDEG 
C 

U  K  I  T  E  (  5 , 1  1  0  )  LUNW  ,  LONE  ,  L  ATVi  , L A  TE ,  LA  TD  ,  LON f) 

110  FUKMATCIOX,  'PARAMETER  CHECK  '  , / ,  y X , 1 7 C 1 H- )  , / / , 

1  bX,'w  LON  =  '  ,  I  b  ,  5X ,  '  u.  LON  =  ',15,/, 

2  5X  ,  '  »  liA'f  =  '  ,  L5,5A,  '  E  LA  f  =  ',15,/, 

3  5X,  'LATK-bA CE=  '  ,  15  ,SX,  'LONE-LONWs  ',1b,/) 

C 

I0IR£C=1 

Ip (  <LA  T.GT.ELAT )  IDlREC=-l 
I0EG1  slUfcG/10 
I  f)  E  G  2  =ID£G-10EG1*10 

c 

UO  200  1  =  1 , LA  TO , 10 EG 
C 

C  CALCULATE  EACH  hIN'S  COORDINATES  wIThIN  GEOGRAPHICALLY 

C  DEFINED  area. 

hLAT  =  LAlWl  +  (.  I0EG*ISIGN*ID1REC*  (  1-1)  ) 

DO  200  0=1 ,  LOUD  ,  1!)EG 
HLU'v  =  LUNHi-t.iDEG*(J-l  J  ) 

C 

C  CALCULATE  CARTESIAN  COORDINATES  OF  VERTICES  FOR  EACH 

C  BIN  AuD  DETERMINE  IE  VERTICES  STRADDLE  PROFILE. 

ICHECK  =0 
pP  =  0 

DO  120  K  =  1 , 2 

BLAll  =HLA  OC  (\-l  )*IUF;G) 

00  120  L=1 , 2 

B L '  "i <“»  1  sdLUi'i  +  (  (  L  -  1  ) *1  DEG ) 

CALu  COORD!  VIJAr,,,fLGN,BLATl,ULUNt,AZMUTH,X,Y) 

C 

C  Y  VALUE  UE  VERT1CF.  *1THIN  ulNDO*? 

IFCAUSCY  )  .  LE . 7 INUUw )  (;U  TU  130 
I CHECKS IChECK+S I GN l 1 . 05 , Y ) 

12<>  CONTINUE 

Ir  (  IAoSUCHECK)  .GE.4)  GO  TO  200 
C 

C  rl  A  w  E  A  iN  APPROPRIATE  dlN  -  SAMPLE  DATA  POINTS. 

130  LAT  =ABS ( ola I ) 

LDim  =  bl,UN 

NEARS=(  LAT/ing;G)  *IOEG 
EARw=(  LOn  /  luEG )  *  1  dEG 
LT 1  = NEARS/ 10 

LT 2  =NEARS-bTi *10 
LNl  =i*EARw/lOO 
LN  2  =in EAR-/10-LN  1*10 
LNi  =NEaR*-(LN1*100+LN2*10) 

C 

ENCODE C 24, 140, LABEL  1 )  H£m(IH£M ) ,LT1 ,LT2,LN1 , LN 2 , LN 3 . IDEG 1 ,1DEG2 
140  F" UK 0 ATC  ’DB2:  1303.3J  B'  ,Al,7ll,  '  .NIS*  ) 

*Kl  t’ECb,  150)  LAttELl 
15u  r'uKMATC'  READING  FILE  -  ',30A1) 

C 

LAdELl C  30 ) =0 

JPE  i  (UNI  T=L«JN2  ,  NABEsLABELl  ,  TY  PE=  '  OLD  '  ,  ACCES5=  '  SEQUENTIAL  '  ,ERR  =  2  00  ) 

1 b 0  REAOC  LUh2 , 17  0 ,  E  N  0  =  1 dO  )  SLaT , S LON , TOPO , F AA 
170  EiJRf- AT  C  2F 10. 4,1  OX  ,E  10.1  ,Eo.l  ,  34X) 

IFCSbUN.LT. 0.0)  SLON=jbO .o+SLON 
C 

C  CALCULATE  CARTESIAN  COORDINATES  OF  DATA  POINTS  IN  BIN. 

CALL  COO HOC  W  LAl’ , KLUN ,SLAT ,SLON , AZHUTH , X , Y ) 

C 

C  ACCEPTABLE  POINT? 

IF C X*CRANGE-X) .LE.O.O)  GO  TU  lbO 

IE CABSCYJ.UT. WINDOW)  GO  TU  160 

C 

C  HAVE  PUIimT  RITHIN  SPECIFIED  WINDOW  -  PLACE  IN  PROJECTION  ARRAYS 

C  IGNORE  REPETITIOUS  POINTS. 

IF C A .E’Q . RPOINT )  GO  TO  lbO 
NPXS  =  H  PT  S  > 1 


A8 


' 


■ 


■ 


JV 

-PLU 

S  V  0  2  ~  5 1  E  1  **  :  2  1  :  SO  2B-JAn-»1 

t  IK,J  .  t  TN 

/  T  RjHLGCi\S/-iiH 

vij<t  y 

■  «P  =:'1P  +1 

U  1  UO 

Pi '5  l  uPTS)=X 

uiLl 

H  P  a  l  N  T  =  X 

GUV  C  VPTS  )  sFAA 

Uioj 

HJP(  *1P  rs  )=TUPO 

''iUH 

IFC.iPTS.EG.20o0)  GO  TO  210 

u  i  u5 

C 

Go  TO  l  h 0 

UiUu 

1  HO 

CLOSE (UN  IT  =  LU w 2 ) 

0  107 

wHITE 15,190)  NP 

Uiub 

iyu 

FuR*ATCiH+,42X# '  POINTS  ACCEPTED  =  ',15) 

u  a  u  y' 

c 

c 

2oO 

CONTINUE 

UilO 

SORT  AuO  OUTPUT  PROJECTION  DATA  FILES. 

c 

c 

210 

CALL  PUbBLE(POS, TOP, GRV , NETS) 

*>  i 

MAKE  POSITIONS  RELATIVE  TO  FEATURE  AT 

C  c 

+-* 

N-  v- 

call  azdis  rc nlat,  nlun  , rl.at, rlon  , differ  ,az) 

DO  220  NC=1,NP1'S 

«UJ 

r* 

22n 

POS(NC)  =PUS(i.C)-OIFEER 

Kil4 

L 

CALL  GCOP  Y  C  POS  ,  TUP  ,  N  PTS  ,  L  A  [*EL2 , 0  ) 

oils 

CALL  GCOP Y l PUS ,GR  V , NPTS , bABEL3  »  0 ) 

V^l  it) 

w  H  l  t  E (  5,230)  u  A  B  E  L  2 , N  p  T  S 

Oil  / 

2  JO 

FORMAT!  1UX TOTAL  NUMBER  OF  PU1NTS  SAMPLED 

1 

1  OX  1  IS  •  ,  Id./) 

WllH 

lFCwPTS.LT. 2D00)  GO  TO  240 

Oil* 

NP  TS  =  0 

GO  TO  ISO 

v>i21 

c 

24  0 

CLOSE  C  L.'<  IT  =  lU'J  2 ) 

ui  22 

L 

90  0 

*t « 1 T  E  (  5 , 4 1 0  ) 

oi2i 

r> 

91o 

FORMAT  I  1  OX  ,  'END  OF  PROCESSING  ') 

oi24 

L 

CALL  EXIT 

ui2b 

Kill) 

PAGE 
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SUBROUTINE  BUBBLE 


AlO 


F  J  K  T  h  A 

0  iV-PMiS  V02-51  r.  21104:00  2 3 -OCT- 

3  u  Li  b  G  h 

.  ir  r  0 

/  1  n  Z  b  GO  C  *  0  /  N  W 

0001 

S  u  0  K  U  U  T  1  0  £  i-}0-Jt'GE(A,vAGl»VAG2,t'O 

c*  + 

BPObGE  3UKT  KUullNfc  lu  MjMEKICAGG 

THE  ASSOCIATED  A  R  K  A  i  6  ,  VAG1,VAG2, 

in  By-'HwrHy  jiih  ahkai  a. 

C** 

0002 

c 

U fbENSIO*  A  c !) , VAGI C 1 ) , VAG2 1 1 ) 

0003 

'i  5  E  G  =0-1 

0004 

DC)  2  0  J  =  1,J 

0005 

f.GC  =0 

0006 

Du  10  I  =  1,.kSEG 

0007 

c 

IF  (  A  l  1  +  1  )  .GT.  All  )  )  GO  TU  10 

00  0  8 

TEWP1  =A (  I J 

00  0  9 

T  E  a  P  2  svALl(I) 

0010 

7EMP3  3VAL2CI) 

00  1  1 

A  (  I  )  =  A  (  I  +  l  ) 

0012 

VAGI  (  U=VAGt  (  Itl  ) 

uu  1  3 

V  A  L  2 ( I)  =  V  A  G  7  C  L +  1  ) 

00  1  4 

A  C  I  4  1  )  =TEMPt 

00  15 

VAGI (  I  +  t  )=  FEMP2 

OOlo 

VAG2C 1+1 )=TEbP3 

0017 

V- 

GdC  =1+1 

00  18 

10 

C •  In  I'INU F 

00  19 

IF ( LOC . GE . 2  )  GO  CO  30 

0O20 

..3EG  =G0C-2 

00  2  1 

20 

O-Oi  ri'Vut 

00  2  2 

L 

30 

R  F  l'  fj  h  0 

002  3 

h  (J 


y  GUKT  AhhAy 
ARE  SOHTfcD 
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I P* 


All 


' 
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SUBROUTINE  COORD 


A12 


21J04I4O 


t'AGt  1 


FURTRrr  I  V  -  P  D  li  3  V02-51F 
COURT) .  Fit*  /  IkiliLiiC'No/  <K 


23-GCl-RO 


OOu  1 


0002 

0  0  0  3 
0004 
0  0  u  5 
0  0  o  o 
000  7 
000  6 

0009 
0010 
00  t  1 
001  2 
0013 

00  1  4 
0015 


3DBRUUTI Nfc  UOOk'H  Hr  La  ,  kFuU,  Kb  A  ,  RLU  ,  THU  TA  ,  XF  ,  YP  ) 

C**  CONVERTS  GEuG*APhIC  C  G  U  K  u  i  n  m  T 1 3  INTO  CAkTESIAI*  CO  UK  DILATES 
C**  *  l  I'H  RESPECT  TO  TuE  PUlul'  kFLA  ,  KFL<U 

C**  AFTER  CONVERSION  ,XP  IS  THE  CGH PO N E N T  IN  DIRECTION  THETA  OKGREES 
C  *  *  Cr.OCRWlSc  From  north,  ip  IS  IN  direction  theta  PLUS  90  DEGREES 

C**  CLOCKWISE. 

C** 

DATA  EARTH, RAO, 1ST ART  /HJ 7 1,0, 57, 29577951 ,0/ 

C 

IF(ISTAKT.EO.l)  GO  TO  10 

CVTDA=EAKTh/K4U 
THETAsTHETA/KAO 
ST  *SI 0 (THETA} 

CT  sCOSC THETA) 

ISTART=1 

C 

10  CVTb0=CVTLA*(C0S(KDA/KAO) ) 
uLA  s(KFLA-RliA)  *CVTLA 
DLO  s-(RFDO-KL.O)  *C  v  TLO 
XP  sDLU*ST-OLA»CT 
YP  sdLA*3T4DDu*CT 

c 

RETURN 

END 
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SUBROUTINE  AZDIST 


A14 
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M4.U1SI . 

t  TU 

/ IK : bLUCKS/  -iK 

uuol 

SUBROUTINE  AZniST  C«liATl ,  Al.OMl  ,ALAT2, ALUN2 . RANGE , THETA  ) 

c** 

c** 

kuUTInE  CALCULATES  THE  RANGE  AND  AZlBUTh 

FROM  (LAT1.L0N1) 

c** 

TO  (UAT2.LON2  )  .  RANGE  IS  IN  KRS.  AZIMUTH 

IN  DEGREES. 

c** 

RANGE  IS  SPHERICAL  DISTANCE  ;  AZIMUTH  IS 

PLANAR  APPROX. 

c** 

c 

DATA  EARTH, R4Q/b37l. 0,5 7. 29577951/ 

UUU  J 

V. 

THETASO.O 

0U(J4 

RanGE=U ,0 

wuob 

R  L  0  N  1  s  A  L  0  N  1 

UOUG 

IFCRLOn l.GT.  180)  RL  UN  1=RUUN  1-360.0 

ooo7 

KLOo  2  = ALUN 2 

OOU8 

r 

1 F ( RLUN2 • GT . I 80 )  RUUN 2=RLUN 2-36 0 . 0 

Oo(j9 

V. 

R  L  A  T 1  s  A  L  A  T 1  /  K  A  D 

UUU) 

KT,Ool  =  RLON  1/KAD 

wn 

KLAi 2=ALAT2/RAU 

uui2 

r* 

RLOn2=KLUN2/H AD 

oi>ii 

V- 

ULUN  = ( R  LUN 1 -RLON  2 ) 

JU  14 

C 

DO  A  T  =lHLAn-KUAi'2) 

wui5 

IF  lOLAT.  nE  .0.0)  GU  TU  20 

UOlo 

KAHGEsEArtTH*A(jSt  oLUN)  *CUS(  C  KLAT  1  +  RLAT2  ) /2 .0) 

ou  1  7 

THETA  =  90 .0 

ouid 

GO  TO  900 

uuiy 

L 

2  o 

ifcduon.ne.o.o)  go  rn  jo 

0  0  20 

k  A  N  G  H  s  A  b  S  (  OLA  C) *EAKTrt 

oo2i 

THETASO . 0 

0022 

IF  (  DU  AT  .GT.  0 . 0  )  THETAs  1  80  .  (j 

ov>23 

r* 

GO  TO  900 

UW24 

L 

3  0 

aLl  SSIN(RUATI) 

0025 

SL2  sSI.N  (  kLiA  T2  ) 

oO  2o 

CL  1  sCUSlRUATl) 

ov>2  7 

C 

CL2  sCOS  (  RU  A  l'2  ) 

o  o  2  b 

COSC  sSLl*SL2+CLl*CL2*CQSC0LUN) 

0U29 

C  saCOSCCOSC) 

Ooju 

R  A  nGE'sABS  (  C*£ARTH ) 

OOJi 

TH£TAsASIN(DLAT/C)*HAD+90.0 

o  o  j  2 

9oo 

IFCDLON .GT.0.0 )  THET As 3b 0 . O-TH ET A 

ooi3 

L 

RETURN 

0OJ4 

EN 1) 

SUBROUTINE  GCOPY 
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FORTR  A 
GCuP  If . 

000  1 


0002 

0003 

0  0  0  4 
0005 
0  0  0  6 
U  00  7 
000  8 
0009 
0010 

0011 
0012 
OOl  J 

0014 
0015 
Oo  l  b 
0017 

()  0  1  8 
0019 
0  0  2  0 
00  21 
0022 

0  0  2  3 
0  0  2  4 
0025 
0  0  2  b 

0  0  2  7 
0028 

0  0  2  9 
0  0  5  0 
0  0  3  1 
0  0  3  2 
0  0  3  3 
0034 

0035 
003b 
0037 
00  3  8 
00  3  9 
o0  40 

0041 
0042 
00  4  3 


.  IV-PLuS  V  u  2  “  5  l  E  2l:u5:09  23-OCT-SO 

•  T  N  /  T  R  !  3  L  U  C  n  S  /  •*  R 


6 U 6 R ul) T I n E  G C 0 P  Y  (  X  ,  G  ,  n  ,  e  lLEfi, IFLAG) 


C** 

C**  HOOriNE  INPUTS  UR  OUTPUTS  DATA  AKhAYS  ACCORDING  to  iflag. 

C**  IFLAG  =  0  -  *OITE  X  AND  G  ARRAYS  INTO  L  A  N>  0  N  T  PRUCT  FUW9AT 

C**  IFLAG  =  l  -  READ  A  AND  G  ARkAYS  INTO  bAMQNT  PRJCT  FORMAT 

C** 

DIMENSION  X ( 1  ) ,G( 1 ) , XX( b) ,  lG(b) 

3 i T  £  FLLEN(l) 

C 

ISTOPsO 

I F ( IFLAG. EG . 1 )  GO  TO  50 
XMNsX(l) 

XHAX=X(N) 

V LAT=0.0 
VLUnsO.O 
A  L M  =  0.0 
C 

OPS’.  (  UrIT=1  ,  0  A4E  =  PILEN  ,  T  Y  P  E  =  ’  N  E '«  '  ,  ACCESS=  '  SEQUENTIAL  '  ) 

•vR  I  TEI  1  ,  10  )  N , VLA  T , VLOn , AAIN , X*AX , AZM 
10  F □RrtAi(I10,5Fl0.J) 

C 

20  L  =0 

ISTARl’  slSTiJP+l 
ISTJP  =  ISTaRT  +  5 
IF  ( 1ST0P.  GT.  .4  3  I3rUP  =  u 
C 

DO  30  JalSTARTrlSTUP 
L  =  L+1 

XX  Cl)  =xiJ) 

IGCL)  =GtJ) 

30  CONTINUE 


C 

**RITE(1,40)  (XX(K)  ,IG(N)  ,K=1, L) 

4  0  FuR*AT(,  IX,  F7. 1  ,Ib,b  7. 1  ,  lb,F7 . 1  ,Ib,F7.  1  ,  Ito,F7 . 1 ,1b  ,F7.1 , 16) 
IF (ISTUP.lt. 4)  GO  TU  20 
GO  TO  500 
C  ’ 

50  OPEN ( UNIT=1 ,  NAMEsEl LEN , TYpE= ' ULD ' ,ACCESS=' SEQUENTIAL'  ) 

READ  (1,10)  '  I  ,  V  L  A  T  ,  V  LUN  ,  AMlN  ,  XMAX,  A2.M 
C 

60  I E  w u  =6 

ISTAR  TsISTOPM 
IS  TOP  =ISTART-t-5 

IFdSTUP.GT.  1)  IEnO=IEND+«-ISTCP 
L  =0 

READ  (1,40)  (XX(I)  ,IG(I)  ,i  =  l,Ifc-fiD) 

C 

DU  70  JsISI'ART,  ISTUP 
L  =L  I 

X(J)  =XX(L) 

<«  (  J  )  =IG(L) 

70  CONTINUE 

IFC ISTCP.LT. N)  GO  Tu  60 

c 

500  CLUSK(UNIT=1 ) 

RETURN 

END 


A 1 7 
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PROGRAM  DISPLY 


AT  8 


FURTP An 

I  V'-PL 

05  V02-51E  2l:Oo:2V  23-UC  T-80  PAGE 

0  I  o  P  U  A  i 

,KTN 

/  r  k  :  h  l  >  c  5  /  8 

0001 

C** 

PROGRAM  DlSPLY 

c** 

ROUTINE  PLOTS  DATA  SIGHT'D  I N  L  A  '-M )  N I  PhOJLCI  FORMAT 

c*  * 

(DEFINED  IN  MRU*  N  buOiv  PnGGRA-S ) . 

c** 

REWRITTEN  FRO*  J.  HOUl.NE'S  P  L  V  E 2  PROGRAM. 

X 0  =  oHIFr  PLuT  XO  INCHES  AuUNG  X  AXIS, 

c** 

Y  0  =  SHIFT  PLOT  xO  INCHES  UP  Y  AXIS 

c** 

XSCALE  =  HUR 1 ZOnTmL  SCALE  (KN/LN) 

c** 

c** 

YSCALE  =  VERTICAL  SCALE  (UnITS/In) 

c** 

c** 

AUT* UR  :  G.D.  KARNtK  (JULY, 1980) 

0002 

DIMENSION  XC  1000)  ,Y(4000 J 

0003 

BYTE  FILEC40) ,ZERU 

0004 

DATA  ZERU/'O'/ 

0005 

u 

10 

»RITE(5,20) 

FORM AT (  '  ENTER  SHIFT  IN  X  AND  Y  DIRECTION  (INCHES)  ',/ 
1  THEN  XSCALe.  AND  YSCALE  ( Jn I TS/ INCH )  '  ) 

0006 

20 

000  7 

ACCEPT  *,  XO.  Y0,  XSC, YSC 
iF( XSC.EU.0.0)  GO  Xu  60o 

0009 

000  9 

xo=-xo*xsc 

0  0  1  0 

YO=-YO*YSC 

0011 

XSCALEsl .O/XSC 

001  2 

c 

:SCALt=l .O/YSC 

001  i 

NR ITE ( 5 , 30 ) 

F  U  K  A  T  (  '  CONTINUOUS  UR  DISCRETE  CURVE  -  LINE  =0,DUT  =1 

0014 

3  0 

0015 

ACCEPT  «,  IDUT 

coin 

*R I  TEC  5, 50) 

0  017 

50 

fohmatc  enter  filename  fur  plotting  ') 

00  1b 

ACCEPT  AC,  Flue 

0  u  1  9 

bO 

FURMAT(40Al) 

0020 

/■» 

FILE(40)sO 

002.1 

CALL  GCOPY (  X  »  Y , NPTS,FILE,+1  ) 

0022 

« 

*RITE(5,o5)  FILE, APTS 

00  2  3 

o  5 

c 

KG R  w  A  T ( 2  X  , 'FILE  -  ',40A1,'  UPEwED. 

1  2X  ,  1 3  ,  '  POINTS  BEING  PLOTTED.  ') 

0  U  2  4 

CALL  SCALF (XSCALE, YSCALE, XO,YO) 

0  0  2  5 

I r  Y  PE  =  " 2 

0  0  2  b 

DU  100  I=1,NPCS 

0027 

XPLUl’xsX  ( I ) 

0026 

YPLUTsY ( I ) 

0029 

IF( 1 DOT, EQ. 0)  GO  TO  70 

0030 

L 

CALL  FCHAR ( XPLUT , YPLOT, 0.001 ,0.001,0.0) 

00  31 

CALL  ANOTAT ( l , 1 H . ) 

00  32 

r* 

GO  TO  100 

00  3  3 

L 

70 

CALL  FPLUTC I  TYPE, XPLUT, (PLOT) 

0034 

I TY  PL  =  0 

U0  35 

100 

r 

CONTINUE 

0  0  3  6 

L 

DU  110  1=1,5 

0037 

YPLuT=(I-3 )/ YSCALE 

0038 

XPL0T=0.1  /XSCALE 

00  39 

CALL  FPL0T(1 ,0.0,  YPLOT) 

0040 

CALL  FPL0T(2«XPLJT, 1PLUT) 

0041 

110 

CON  rif.UE 

0042 

L 

XPLGT=-0.2/XSCALE 

0043 

CALL  FCHAR (XPLUT, 0.0, 0.1, 0.1, 0.0) 

0044 

/•* 

CALL  A.NOTATC  1  ,Z£RO  ) 

0045 

L 

YPLQT= 1.0/ YSCALE 

0O4o 

CALL  FCHAR (0.0, YPLOT, 0.1, 0.1, 0.0) 

0047 

r* 

CALL  AnQTAT(40,FILE) 

0048 

L 

XPL0T=1.5/XSCALE 

0049 

Y P  LOT= 1 . 0 /YSCALE 

0050 

CALL  FCHAR (XPLUT, YPLOT, 0.1,0. 1 ,0.0) 

A19 
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DISPLm t , 

IV-PLUS  M2-51E  21:0t>:  27  23-UCT-hu  PAGE  2 

,  t'TN  /'I  a  :  eLuCr^S/  *H 

005  1 

0052 

0<>5.H 

EoCUDEChG, 120, FILE)  XSC,ii>C 

120  KOHr-AT  (yHXSCALE  =  ,  F 1  0  .  i  ,  2X  ,  9H  YSCALE  =  ,F10.1) 
call.  ALOTATC 40 , FILE ) 

C 

0054 

0055 

CALL  FPLUTC 1 ,0.0 ,0.0) 

GO  TU  10 

oo 

oo 

crc* 

'-JO 

C 

500  CALL  EXIT 

F.  O  O 

t 
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Appendix  B 


Printouts  from  programs  BINNER  and  PBGT 
showing  examples  of  their  application  . 


RUN  BINNER 
00:00:17 

ENTER  DATA  AREA  -  WLON , ELON ,SLAT ,  NLAT 

-002, +002, -01 ,+01 

N  x  N  BLOCK  SIZE  (INTEGER  DEGREES) 

01 

CREATE  FILES  FOR  1  x  1  BLOCKS?  YES=1 
1 

PARAMETER  CHECK 


W  LON  =  358  E  LON  =2 
S  LAT  1  N  LAT  =  1 
LATN-LATS=  3  LONE-LONW  =  5 


(  1  ) BS01 35801 . N I B ;  1 

(  2  ) BS01 35901 .NIB ;1 

(  3  ) BS01 00001 . N I B ;  1 

(  4  ) BS01 001 01 . NIB ;1 

(  5  ) BS01 00201 . N I B ; 1 

(  6  ) BN003 5801 . N I B ; 1 

(  7  ) BN0035901 .NIB ;1 

(  8  )BN0000001 . N I B ; 1 

(  9  ) BN00001 01 . N I B ; 1 

(  10  ) BN0000201 . N I B ; 1 

(  11  ) BN01 35801 .NIB ;1 

(  12  )BN01 35901 . N I B ; 1 

(  13  ) BN01 001 01 . NIB ;1 

(  14  ) BN01 00201 . N I B ; 1 

TYPE  RAW  DATA  FILE  NAME  -  30A1 
APPALACHIAN.DAT 


B1 


& 

' 


RUN  BINNER 

-001.1,001.1,-1.1,1.1 

ENTER  DATA  AREA  -  WLON , ELON , SLAT , NLAT 

N  x  N  BLOCK  SIZE  (INTEGER  DEGREES) 

1 

CREATE  FILES  FOR  1  x  1  BLOCKS?  YES=1 
1 

PARAMETER  CHECK 


w 

LON  =  358 

E  LON 

s 

LAT  =  2 

N  LAT 

LATN-LATS=  5 

LONE-LONW 

( 

1 

) BS02  35801 . NIB ;1 

( 

2 

) BS0235901 . NIB ;  1 

( 

3 

)BS0200001 . N I B ; 1 

( 

4 

) BS02001 01 . N I B ;  1 

( 

5 

)BS0200201 . N I B ; 1 

( 

6 

) BS01 35801 . N I B ; 1 

( 

7 

) BS01 35901 . N I B ; 1 

( 

8 

) BS01 00001 . N I B ; 1 

( 

9 

) BS01 001 01 . N I B ;  1 

( 

10 

) BS01 00201 . N I B ; 1 

( 

11 

)BN0035801 . NIB ;1 

( 

12 

) BN003  5901 . N I B ;  1 

( 

13 

)BN0000001 . N I B ; 1 

(. 

14 

)BN00001 01 . N I B ;  1 

( 

15 

)BN0000201 . N I B ;  1 

( 

16 

) BN01 35801 . N I B ; 1 

( 

17 

) B  N  0 1 35901 . N I B ;  1 

( 

18 

) BN01 00001 .NIB  ;1 

B2 


RUN  BINNER' 

23:54:50 

ENTER  DATA  AREA  -  WLON , ELON , SLAT , NLAT 
-007, +012, -03, +18 
N  x  N  BLOCK  SIZE  (INTEGER  DEGREES) 

10 

CREATE  FILES  FOR  10  x  10  BLOCKS?  YES=1 
1 

PARAMETER  CHECK 


W  LON 

350 

E  LON 

12 

S  LAT 

10 

N  LAT 

28 

LATN-LATS= 

48 

LONE-LONW= 

32 

(  1  )BS1035010.NIB;1 

(  2  ) BS 1 00001 0 . N I B ; 1 

(  3  )BS1001010.NIB;1 

(  4  )BS1002010.NIB;1 

(  5  ) BN003501 0 . NIB ; 1 

(  6  )BN000001 0 .NIB ;1 

(  7  )BN0001010.NIB;1 

(  8  )BN0002010.NIB;1 

(  9  )BN1035010.NIB;1 

(  10  ) BN1 00001 0 . N I B ; 1 

(  11  ) BN1 001 01 0 .NIB ;1 

(  12  ) BN 1 00201 0 . N I B ; 1 

(  13  )BN2035010.NIB;1 

(  14  )BN2000010.NIB;1 

(  15  ) BN2001 010.NIB;! 

(  16  )BN2002010.NIB ;1 

(  17  )BN3035010.NIB;1 

(  18  )BN3000010.NIB;1 

(  19  ) BN3001 010.NIB;1 

(  20  )BN3002010.NIB;1 

TYPE  RAW  DATA  FILE  NAME  -  30A1 


. 

. 

■ 


RUN  BINNER 
20:53:59 

ENTER  DATA  AREA  -  WLON , ELON , SLAT , NLAT 
-125,-118,48,49 

N  x  N  BLOCK  SIZE  (INTEGER  DEGREES) 

1 

CREATE  FILES  FOR  1  x  1  BLOCKS?  YES=1 
1 

PARAMETER  CHECK 


W  LON  =  235  E  LON  =  242 
S  LAT  =  48  N  LAT  =  49 
LATN-LATS=  2  LONE-LONW=  8 


(  1  ) BN482  350 1 . N I B ; 1 

(  2  )BN4823601 . N I B ; 1 

(  3  ) BN482 3701 .NIB ;1 

(  4  )BN4823801 . N I B ; 1 

(  5  ) BN482 3901 . N I B ; 1 

(  6  )BN4824001 .NIB; 1 

(  7  ) BN48241 01 . N I B ; 1 

(  8  ) BN4824201 . N I B ; 1 

(  9  ) BN4923501 . N I B ; 1 

(  10  ) BN4923601 . N I B ; 1 

(  11  ) BN492  3  701 . N I B ; 1 

(  12  ) BN4923801 . N I B ; 1 

(  13  ) BN4923901 . N I B ; 1 

(  14  ) BN4924001 . N I B ; 1 

(  15  ) BN49241 01 . NIB ;1 

(  16  ) BN4924201 . N I B ; 1 

TYPE  RAW  DATA  FILE  NAME  -  30A1 
DB2:USHRT. FAE 


B4 


, . 


DIR  * . N I B ; * 


DIRECTORY  DB1: [303,26] 
27-MAR-80  21:31 


BN4823501 . N I B ; 1 

7 

27-MAR-80  20:54 

BN4823601 .NIB;1 

11 

27-MAR-80  20:54 

BN4823701 . N I B ; 1 

51 

27-MAR-80  20:54 

BN4823801 . NIB ;1 

25 

27-MAR-80  20:54 

BN4823901 .NIB;1 

42 

27 -MAR-80  20:54 

BN4824001 .NIB;1 

53 

27-MAR-80  20:54 

BN4824101 . N I B ; 1 

53 

27-MAR-80  20:54 

BN4824201 . N I B ; 1 

1 

27-MAR-80  20:54 

BN4923501 . N I B ; 1 

1 

27-MAR-80  20:54 

BN4923601 . N I B ; 1 

1 

27-MAR-80  20:54 

BN4923701 .NIB;1 

1 

27-MAR-80  20:54 

BN4923801 .NIB;1 

0 

27-MAR-80  20:54 

BN4923901 . N I B ; 1 

0 

27-MAR-80  20:54 

BN4924001 . N I B ; 1 

0 

27-MAR-80  20:54 

BN49241 01 .NIB;1 

1 

27-MAR-80  20:54 

BN4924201 . N IB ;  1 

0 

27-MAR-80  20:54 

TOTAL  OF  247/285  BLOCKS  IN  16  FILES 
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Figure  2.  An  example  of  the  use  of  program  DISPLY  using  as  input,  the  profile 
constructed  from  figure  1  using  program  PBGT.  These  files  are 
stored  as  sequential,  PRJCT  format  files  (as  defined  by  BROWN 
BOOK  program  PRJCT). 
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